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Abstract 



This report describes the first search for top squark pair production in the 
channel titi bbxtxi ~^ ee+jets+^y using 74.9 ± 8.9 ph-^ of data col- 
lected using the D0 detector. A 95% confidence level upper limit on a ■ B is 
presented. The limit is above the theoretical expectation for cr • B for this pro- 
cess, but does show the sensitivity of the current D0 data set to a particular 
topology for new physics. 

PACS numbers: 14.80.Ly, 13.85.Rm 
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Supersymmetry (SUSY) is a fundamental space-time symmetry relating bosons and 
fermions Supersymmetric extensions to the Standard Model (SM) feature undiscov- 
ered superpartners for every SM particle — for example, there is a scalar quark (squark) 
for each of the two degrees of freedom for the spin 1/2 quarks. In most SUSY models, the 
masses of the squarks are approximately degenerate except for those of the top squarks. 
Due to large top family Yukawa interactions, the lighter top squark mass eigenstate (ti) can 
have a much lower mass than the other squarks Top squarks will be pair produced at 
the Tevatron; each will then decay into the lightest chargino Xi ^ind a b quark if that decay 
is kinematically allowed. If m~± is greater than the mass of the top squark, the decay will 
proceed through slepton channels. If the sleptons are also heavier than the top squark, it 
will decay to a charm quark and the lightest supersymmetric particle (LSP or Xi) with a 
100% branching fraction. 

The Xi decays to li^iXi Qq'Xi- Under i?-parity conservation, the lightest neutralino 
is assumed to be stable and escapes detection, resulting in missing transverse energy 
Thus, top squarks, pair-produced at the Tevatron, result in final states similar to those of 
top quarks. However, as the decay of the chargino is to three particles, the decay products 
tend to be softer than those of the W boson. Since the xf decay is almost always dominated 
by virtual W exchange, the branching fractions are expected to be be very close to W boson 
leptonic and hadronic decay branching fractions 0. 

The results of a search for ti cXi have been published by the D0 Collaboration [Q. 
Model independent lower limits on the masses of the top squark and lightest chargino have 
been set using the measured width of the Z boson and are approximately 45 GeV/c^ [Q. 
Within the framework of the Minimal Supersymmetric Standard Model (MSSM) the 
current limits on the pair production of charginos (which depend on the assumed value 
of the common scalar mass mo) from LEP at ^/s = 130, 136, and 161 GeV 0,0], lead to 
m~± > 62.0 - 78.5 GeV/c^ at the 95% CL 0. The analysis described below is independent 
of the MSSM and supergravity ^ frameworks, instead depending only on the masses of the 
top squark, the lightest chargino, and the lightest neutralino, and on the branching fractions 
of the chargino decay. 

Previous phenomological studies have considered final states with a single lepton -|- jets 
+ and two leptons + jets + J^^ [@]- These studies, which used ISAJET events smeared 
by typical detector resolutions, indicated that the single lepton channel cannot be studied at 
the Tevatron without excellent 6-tagging capability due to the enormous background from 
W boson production. However, they did indicate that analysis of the dilepton channels (ee, 
efi and /i/x) could lead to a limit on the mass of (or discovery of) the top squark at the 
Tevatron using the current data set. 

This report describes the first search for the decay ti — > bxf in the channel tiii — ^ ee+ 
jets + ]^rp using 74.9 ± 8.9 pb~^ of data. The data were collected at the Fermilab Tevatron 
at y/s = 1.8 TeV during 1994-1995. The D0 detector and data collection system are 
described in detail in Ref. |10]. The detector consisted of three major subsystems: a uranium- 
liquid argon calorimeter, central tracking detectors (with no central magnetic field), and a 
muon spectrometer. Electrons were identified by their longitudinal and transverse shower 
profiles in the calorimeter and were required to have a matching track in the central tracking 
chambers. In this analysis, they were restricted to have pseudorapidity |?7| < 2.5 and to be 
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isolated from other energy depositions in the event. Jets were reconstructed using a cone 
algorithm of radius IZ = ^ (A</))2 + (Arj)'^ = 0.7 with \r]\ < 4.0. The $]rp was determined 
from energy deposition in the calorimeter for \t]\ < 4.5. 

The acceptance for top squark events was calculated for a range of top squark and 
chargino masses using the ISAJET event generator and a detector simulation based on the 



GEANT program [^. Samples were generated with top squark masses between 55 and 75 
GeV/c^ with m~± between 47 and 68 GeV/c^, depending on nij^. The mass of the lightest 
neutralino was set to the supergravity-motivated value 

The signature for titi ^XiXi is two electrons, one or more jets, and Kinematic 
distributions for {mj^,m~±) = (65,47) GeV/c^ are shown in Fig. ^ This analysis was 
restricted to events selected using a trigger which required one electromagetic cluster with 
transverse energy > 15 GeV, one jet with Ej, > 10 GeV, and > 14 GeV. Other 
kinematic quantities used to discriminate against background are the invariant mass of the 
two electrons rriee and E^'^ = E^ + E^ + (defined in Ref. as bigness). 

Cut optimization was done using the RGSEARCH program. RGSEARCH uses a 
modified grid search based on Monte Carlo (MC) signal events and background samples to 
optimize event selection. In this study, the MC signal samples described above and the MC 
physics background samples listed in Table | were used. Several combinations of selection 
criteria were explored starting with the thresholds imposed by the trigger conditions. The 
final selection criteria are summarized in Table ||. Other combinations included require- 
ments on the Et of a second jet and/or the azimuthal angle between the two electrons. 
These combinations increased the signal to background ratio, but reduced the signal effi- 
ciency significantly. Values for the upper limits on rriee and i?™™ were fixed while running 
RGSEARCH. The cut on rriee was used to remove Z —>■ ee events and that on i?^"™ to remove 
ti ee-t-jets-t-^-j^ events. Distributions of E^^ for top squark production with (m^^,m~±) 
= (65,47) GeV/c^ and Monte Carlo top quark production are shown in Fig. 0. 

Signal detection efficiency was restricted by the reconstruction and identification of low 
Et electrons. Only approximately 15% of Monte Carlo events with (mj^,m~±) = (65,47) 
GeV/c^ had two reconstructed electromagetic clusters (with an associated track) with E^ > 
16 GeV and E^ > 8 GeV. In addition, the identification efficiency for two electrons, one 
with Et = 8 GeV and one with Et = 16 GeV, was approximately 40%. It is, however, 
essential to include the second electron in the selection criteria to avoid being overwhelmed 
by W boson events. 

Physics backgrounds were estimated by Monte Carlo simulation or from a combination of 
Monte Carlo and data. The instrumental background from jets misidentified as electrons was 
estimated entirely from data |jl3[| using the jet misidentification probability for the electron 
identification and kinematic cuts used in this analysis ((6.5 ± 1.3) x 10~^). Four physics 
backgrounds were considered in this study: tt production with a top quark mass of 170 
GeV/c^, WW production, and Z boson production with final states resulting in dielectrons. 
The contribution to the background from individual channels is given in Table |l|. The total 
predicted background is 4.4 ± 0.8 events. 

After application of the cuts to the data sample, two events remained. Given no observed 
excess of events above the expected background, we set a 95% CL upper limit on a ■ B 
using a Bayesian approach with a fiat prior distribution for the signal cross section. 
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The statistical and systematic uncertainties on the efficiency, the integrated luminosity, 
and the background estimation were included in the limit calculation with Gaussian prior 
distributions. The resulting upper limit on cr • i? as a function of nij^ with fixed m~± = 47 
GeV/c^ is shown in Fig. ^ along with the predicted a ■ B. The choice of m~± = 47 GeV/c^ 
allows the widest range of m-^^. As can be seen, no limit on mj^ can be set. The situation is 
similar for increased m~±. 

Although the recent results on chargino pair production from LEP limit the likelihood for 
a light top squark to decay to a quark and a chargino within the MSSM, the a-B limit curve 
shown in Fig. ^ indicates the level of sensitivity in the current D0 data set to a particular 
topology for new physics: pair production of new particles which decay into leptons, jets, and 
non-interacting particles. Such a new particle, with a top-like signature, could be detectable 
in the current data set down to a production cross section times branching ratio of order 10 
pb. 
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and acknowledge support from the Department of Energy and National Science Foundation 
(U.S.A.), Commissariat a L'Energie Atomique (France), Ministries for Atomic Energy and 
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Science and Education (India), Colciencias (Colombia), CONACyT (Mexico), Ministry of 
Education and KOSEF (Korea), CONICET and UBACyT (Argentina), and the A. P. Sloan 
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TABLES 



TABLE L Background contributions from individual channels. 



Background Channel 


Number of Events 


tt(170) ^ ee 


0.03 ± 0.01 


WW^ee 


0.02 ± 0.01 


Z ^ee 


0.09 ±0.01 


Z ^ TT ^ ee 


0.67 ±0.13 


Misidentification 


3.6 ±0.8 


Total 


4.4 ±0.8 



TABLE II. Kinematic cuts. E^™ is defined in the text. 

E!j.^ > 16 GeV 
Elf > 8 GeV 
E^j^ > 30 GeV 
> 22 GeV 
ruee < 60 GeV/c2 
E^"^ < 90 GeV 
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FIGURES 
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FIG. 1. Kinematic distributions for (mr ,m~±] 



(65,47) GeV/c^ 
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FIG. 2. Distributions of E?^'^ for (a) {mj^ ,m~± ) = (65,47) GeV/c^ and (b) top quark production 
with mt = 170 GeV/c^. 



9 




55 60 65 70 75 

(GeV/c^) 

FIG. 3. Our 95% confidence level upper limit on a ■ B as a function of mr for m~± = 47 
GeV/c^. Also shown are the predicted values from IS A JET. 



10 



